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Abstract 

We measure the spectral dimension of universes emerging from nonperturbative 
quantum gravity, defined through state sums of causal triangulated geometries. 
While four-dimensional on large scales, the quantum universe appears two-di- 
mensional at short distances. We conclude that quantum gravity may be "self- 
renormalizing" at the Planck scale, by virtue of a mechanism of dynamical dimen- 
sional reduction. 
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Quantum gravity as an ultraviolet regulator? 

A shared hope of researchers in otherwise disparate approaches to quantum grav- 
ity is that the microstructure of space and time may provide a physical regulator 
for the ultraviolet infinities encountered in perturbative quantum field theory. 
The outstanding challenge is to construct a consistent quantum description of 
this highly nonperturbative gravitational regime that stands a chance of being 
physically correct. 

Slow progress in the quest for quantum gravity has not hindered speculation 
on what kind of mechanism may be responsible for resolving the short- distance 
singularities. A recurrent idea is the existence of a minimal length scale, often 
in terms of a characteristic Planck-scale unit of length in scenarios where the 
spacetime at short distances is fundamentally discrete. An example is that of 
so-called loop quantum gravity, where the discrete spectra of geometric operators 
measuring areas and volumes on a kinematical Hilbert space are often taken 
as evidence for fundamental discreteness in nature 13 1 Other quantization 
programs for gravity, such as the ambitious causal set approach [5], postulate 
fundamental discreteness at the outset. 

The alternative we will advance here is based on new results from an anal- 
ysis of the properties of quantum universes generated in the nonperturbative 
and background-independent CDT (causal dynamical triangulations) approach 
to quantum gravity. As shown in E] , they have a number of appealing macro- 
scopic properties: firstly, their scaling behaviour as function of the spacetime 
volume is that of genuine isotropic and homogeneous four-dimensional worlds. 
Secondly, after integrating out all dynamical variables but the scale factor a(r) 
in the full quantum theory, the correlation function between scale factors at dif- 
ferent (proper) times r is described by the simplest minisuperspace model used 
in quantum cosmology. 

We have recently begun an analysis of the microscopic properties of these 
quantum spacetimes. As in previous work, their geometry can be probed in a 
rather direct manner through Monte Carlo simulations and measurements. At 
small scales, it exhibits neither fundamental discreteness nor indication of a min- 
imal length scale. Instead, we have found evidence of a fractal structure (see [7] , 
which also contains a detailed technical account of the numerical set-up). What 
we report on in this letter is a most remarkable finding concerning the universes' 
spectral dimension, a diffeomorphism-invariant quantity obtained from studying 

1 Although this may be an attractive scenario, it is quite possible that the discreteness (and 
associated finiteness) results of loop quantum gravity will not survive the transition to the 
physical sector of the theory, which also solves the Hamiltonian constraint. One may also 
question the genericity of the "emergence" of such discreteness, which is a consequence [3] of 
the insistence that one-dimensional, unsmeared Wilson-line variables be represented by finite 
operators in the quantum theory. 
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diffusion on the quantum ensemble of geometries. On large scales and within 
measuring accuracy, it is equal to four, in agreement with earlier measurements 
of the large-scale dimensionality based on the scale factor. Surprisingly, the spec- 
tral dimension turns out to be scale-dependent and decreases smoothly from four 
to a value of around two as the quantum spacetime is probed at ever smaller dis- 
tances. This suggests a picture of physics at the Planck scale which is radically 
different from frequently invoked scenarios of fundamental discreteness: through 
the dynamical generation of a scale-dependent dimensionality, nonperturbative 
quantum gravity provides an effective ultraviolet cut-off through dynamical di- 
mensional reduction. 



The spectral dimension 

A diffusion process on a <i-dimensional Euclidean geometry with a fixed, smooth 
metric g a b(0 is governed by the diffusion equation 

■j^K g (Z,to*) = A B K g (Z,fa<T), (1) 

where a is a fictitious diffusion time, A g the Laplace operator corresponding to 
9ab(0 and K g (£, £ ; cr) the probability density of diffusion from £ to £ m diffusion 
time a. We will consider processes which are initially peaked at some point £o> 

K g (Z,faa = 0) = ?}LJ& . (2) 

A quantity that is easier to measure than K g in numerical simulations is the 

average return probability 

P g (a) := iy d d {y/tetgj£)K g (t,t-a), (3) 

where V = J d d ^ ydet g(£) is the spacetime volume. 

For an infinite flat space, the solution to eq. is simply given by 

e -4(£>?o)/4<r 

Kg(£,fa<r) = ( 47nT W2 ' 9ab(0=$ab, (4) 

where c? 9 (£,£o) denotes the geodesic distance between £ and £o- It follows that 
yfa is an effective measure of the linear spread of the Gaussian at diffusion time 
a. Because of P g (o~) = l/a d ^ 2 in the flat case, we can extract the dimension d of 
the manifold by taking the logarithmic derivative, 

d log o 
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independent of a. 

For curved spaces and/or finite volume V one can still use eq. (jSJ) to extract the 
dimension, but there will be corrections for sufficiently large a. For finite volume 
in particular, P g (o~) goes to one for a ^> V 2 l d since the zero mode of the Laplacian 
— A g will dominate the diffusion in this region. For a given diffusion time a the 
behaviour of P g {o) is determined by eigenvalues A n of — A g with A n < 1/cr, and 
the contribution from higher eigenvalues is exponentially suppressed. Like in the 
flat case, large a is related to diffusion which probes spacetime at large scales, 
whereas small a probes short distances. 

We will use the return probability to determine an effective dimensionality 
of quantum spacetime, which in our nonperturbative path integral formulation 
amounts to studying diffusion on an entire ensemble of curved, Euclidean(ized) ge- 
ometries. Since the return probability P g {cr) in (J3J) is invariant under reparametri- 
zations, one can define the quantum average -Py(cr) of P g (cr) over all equivalence 
classes [g a b] of metrics with a given spacetime volume V by 

Pv(ct) = ±- jv[g ab ] e - s ^5( J d 4 £^tg~ - V) P», (6) 

where SE(g a b) is the (Euclidean) Einstein-Hilbert action of g a b(0 and Zy the 
partition function (path integral) for geometries with fixed volume V. 

It is straightforward to generalize the diffusion process and its associated 
probability density to the piecewise linear geometries that appear in the path 
integral of CDT [7J- In analogy with the ordinary path integral for a particle, 
one would expect that in the continuum limit a typical geometry in the quantum 
ensemble is still continuous, but nowhere differentiable; for example, it could be 
fractal. In fact, diffusion on fractal structures is well studied in statistical physics 
jH], and there the return probability takes the form 

P„M = ^ F (_^_), (7 ) 

where N is the "volume" associated with the fractal structure and D$ is the 
so-called spectral dimension, which is not necessarily an integer. An example 
of fractal structures are branched polymers, which generically have a spectral 
dimension D s =A/3 EDI- 

For the nonperturbative gravitational path integral defined by CDT, we expect 
the same functional form ((Zj), where iV now stands for the discrete four- volume, 
given in terms of the number of four-simplices of a triangulation. As above in 
(JHJ), we can now extract the value of the fractal dimension D$ by measuring the 
logarithmic derivative, 

d\ogP N (a) 

P>s{or) = -2 — , (8) 

d log a 
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as long as the diffusion time a is not much larger than N 2 ^ Ds . From previous 
numerical simulations of 2d quantum gravity in terms of dynamical triangulations 
P"H IT3| we already know that eq. (JSJ) is not reliable for arbitrary diffusion 
times. The observed behaviour of Ds{cr) for a given triangulation will typically 
exhibit irregularities for the smallest a, caused by the lattice discretization 2 , and 
then enter a long and stable regime where the spectral dimension is independent 
of a, before finite-size effects start to dominate and Ds(cr) goes to zero. 

Measuring the spectral dimension 

In the CDT approach, quantum gravity is defined as the continuum limit of a 
regularized version of the nonperturbative gravitational path integral [TH ITo] . 
The set of spacetime geometries to be summed over is represented by a class of 
causal four-dimensional piecewise flat manifolds ("triangulations"). Every mem- 
ber T of the ensemble of simplicial spacetimes can be wick-rotated to a unique 
Euclidean piecewise flat geometry, whereupon the path integral assumes the form 
of a partition function 



where Ct is a combinatorial symmetry factor and Se{T) the Euclidean Einstein- 
Regge action of the triangulation T ^3]. All geometries share a global, discrete 
version of proper time. In the continuum limit, the CDT time r becomes propor- 
tional to the cosmological proper time of a conventional minisuperspace model 



We extract information about the continuum limit of the theory by Monte 
Carlo simulations and a finite-size scaling analysis of (JHJ). Further details about 
the updating moves for the geometry and the numerical set-up can be found in ITol 
and [Zj respectively. For computer-technical reasons we keep the total spacetime 
volume N approximately fixed during simulations. All the results obtained can 
be related by an inverse Laplace transform to those for the geometric ensemble 
with the volume constraint absent. 

Our measurements of the spectral dimension were performed in the phase of 
the statistical model which generates a quantum geometry extended in both 

2 Often the behaviour of Pjsr(cr) for odd and even numbers of diffusion steps a will be quite 
different for small a and merge only for a sa 20 — 30. - The origin of this asymmetry is illustrated 
by the (extreme) example of diffusion on a one-dimensional piecewise straight space, where the 
return probability simply vanishes for any odd number <r of steps. 

3 The proper time r of CDT - although invariantly defined from a geometric point of view 
- does not necessarily coincide with a physical time appearing explicitly in the formulation of 
physical observables, for example, two-point functions. 




(9) 
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space and time, with large-scale four- dimensional scaling properties 0I7]. As dis- 
cussed in 0, the quantum universe we generate is a "bounce" in the Euclidean 
sector of the theory. It has a characteristic shape when we plot its three-volume 
(equivalently, its scale factor a(r)) as a function of proper time r. The universe 
starts out with a minimal three-volume, increases to a maximum, and then de- 
creases in a symmetric fashion back to a minimal value. As shown in the 
underlying dynamics of the scale factor a(r) solves the equation of motion of a 
simple minisuperspace action for a(r), and is that of a bounce, with total Eu- 
clidean spacetime volume determined by the volume at which the simulation is 
performed. It is on an ensemble of geometries of this type that we made our 
measurements. Since we are interested in the bulk properties of the diffusion, we 
always started the process from a simplex in the constant-time slice where a(r) is 
maximal. The simulations to determine the universe's spacetime spectral dimen- 
sion were performed for geometries of proper-time extension t = 80 and a discrete 
volume of up to approximately A/" = 181.000 four-simplices, and the diffusion was 
followed for up to cr m ax = 400 time steps. 

Fig. summarizes our measurements of the spectral dimension Ds{<j) at the 
maximal spacetime volume, extracted as the logarithmic derivative (jSJ) from a 
discrete implementation of the diffusion process. The (envelopes of the) error 
bars represent the errors coming from averaging over 400 different measurements 
of diffusion processes, performed for independent starting points and statistically 
independent configurations T generated by the Monte Carlo simulation. As ob- 
served previously in other systems of random geometry, we have found a different 
behaviour of Ds(<j) for odd and even (discrete) diffusion times a for small a (cf. 
footnote|2]). In order to eliminate this short- distance lattice artifact, we have only 
included the region a > 40 for which the odd and even curves coincide, both in 
Fig. ^ and in determining the spectral dimension. 

The data points along the central curve in Fig.Qrepresent our best approxima- 
tion to Ds(cr) in the limit of infinite spacetime volume. Their monotonic increase 
as function of a indicates that we have not yet reached the region where finite- 
volume effects dominate (in the form of the constant mode of the Laplacian). 
The remarkable feature of the behaviour of the spectral dimension illustrated in 
Fig. n is that it is qualitatively different from what has been observed in similar 
systems up to now, be it in two-dimensional Euclidean quantum gravity with or 
without matter ^T] or for the spatial hypermanifolds of our present CDT set- 
up in four dimensions |7j. In these cases, immediately following the region of 
even-odd asymmetry for small a, Ds(<j) stabilizes in a horizontal line, indicat- 
ing the presence of a single value Dg characterizing the spectral dimension of 
the system, independent of the scale at which the diffusion process probes the 
geometry. 

Apparently, this is not the case for the spectral spacetime dimension in quan- 
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Figure 1: The data points along the central curve show the spectral dimension Dg(a) 
of the universe as function of the diffusion time a. Superimposed is a best fit, the 
continuous curve Ds(cr) = 4.02 — 119/(54+cr). The two outer curves quantify the error 
bars, which increase linearly with a, due to (JHJ). (Measurements taken for a quantum 
universe with 181.000 four-simplices.) 



turn gravity defined by CDT. The measurements shown in Fig. ^ indicate that 
Ds(<j) changes with the scale probed. In order to quantify this scale dependence 
we have attempted a variety of fits in the available data range a G [40,400]. 
Among curves with three free parameters of the form consi.+asymptotic form 4 , 
a fit of the form 

_ 2 i}°im =a _j L _ (10) 

a log a a + c 

agrees best with the data. In Fig. [TJ the curve 

D s (a) = 4.02 - (11) 
54 + a 

has been superimposed on the data, where the three constants were determined 
from the entire data range a G [40,400]. Although both b and c individually 
are slightly altered when one varies the range of a, their ratio b/c as well as the 
constant a remain fairly stable. Integrating relation ()10p. we have 

P{*) - 3T7^ (12) 



(W^l + c/a) 6 / 



4 The two alternatives considered were a— be ca and a~b/a c 
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implying a behaviour 




a 



a/2 



for large cr, 



P(a) 



(13) 



(a-6/c)/2 



for small cr. 



Our interpretation of eqs. (fT2J) and (fT3j) is that t/ie quantum geometry generated 
by CDT does not have a self-similar structure at all distances, but instead has a 
scale- dependent spectral dimension which increases continuously from a—b/c to a 
with increasing distance. 

Taking into account the variation of a in eq. (fTUjl when using various cuts 
[cminj °max] f° r the range of cr, as well as different weightings of the errors, we 
obtain the asymptotic value 



which means that the spectral dimension extracted from the large-cr behaviour 
(which probes the long-distance structure of spacetime) is compatible with four. 
On the other hand, the "short-distance spectral dimension" , obtained by extrap- 
olating eq. (|12|) to cr — > is given by 



and thus is compatible with the integer value two. 
Discussion 

The continuous change of spectral dimension described in this letter constitutes 
to our knowledge the first dynamical derivation of a scale-dependent dimension in 
full quantum gravity. 5 It is natural to conjecture it will provide an effective short- 
distance cut-off by which the nonperturbative formulation of quantum gravity 
employed here, causal dynamical triangulations, evades the ultraviolet infinities 
of perturbative quantum gravity. Contrary to current folklore (see ^7] for a 
review), this is done without appealing to short-scale discreteness or abandoning 
geometric concepts altogether. 

Translating our lattice results to a continuum notation requires a "dimensional 
transmutation" to dimensionful quantities, in accordance with the renormaliza- 
tion of the lattice theory. Because of the perturbative nonrenormalizability of 
gravity, this is expected to be quite subtle. CDT provides a concrete frame- 
work for addressing this issue and we will return to it elsewhere. However, since 

5 In the so-called exact renormalization group approach to Euclidean quantum gravity, a 
similar reduction has been observed recently in an Einstein-Hilbert truncation |16| . 




(14) 




(15) 
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a from (|TJ can be assigned the length dimension two, and since we expect the 
short-distance behaviour of the theory to be governed by the continuum gravita- 
tional coupling G^, it is tempting to write the continuum version of (JTUJ) as 

Pv{<?) ~ -r z t „ i , (16) 

K ' a 2 1 + const. ■ G N /a v ' 

where const, is a constant of order one. The relation (fTBj) describes a universe 
whose spectral dimension is four on scales large compared to the Planck scale. 
Below this scale, the quantum-gravitational excitations of geometry lead to a 
nonperturbative dynamical dimensional reduction to two, a dimensionality where 
gravity is known to be renormalizable. 
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